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Different choices exist for the renormalisation group resummation in the determination of a s from hadronic r 
decays: namely fixed-order (FOPT) and contour-improved perturbation theory (CIPT). The two approaches lead 
to systematic differences in the resulting a s . On the basis of a model for higher-order terms in the perturbative 
series, which incorporates well-known structure from renormalons, it is found that while CIPT is unable to account 
for the fully resummed series, FOPT smoothly approaches the Borel sum. Employing the model to determine a„ 
yields a a (M T ) = 0.316 ± 0.006, which after evolution leads to a s {Mz) = 0.1180 ± 0.0008. 



1. INTRODUCTION 

Due to its particular mass of M T « 1.8 GeV, 
the t lepton constitutes an excellent system to 
study QCD at low energies, as in about 65% of 
all cases it decays into hadrons, while the QCD 
description remains largely perturbative. In the 
seminal work [1 , the framework for a precise de- 
termination of the QCD coupling a s from the to- 
tal r hadronic width 



R T 



T\t- 



hadrons f T (7)] 



r[r- 



3.640(10), (1) 



was developed, while afterwards invariant mass 
distributions were incorporated into the analy- 
sis as well |2l3j . The most recent study of the 
ALEPH spectral function data on the basis of 
the final full LEP data set yielded a s (M T ) — 
0.344 ± 0.005 exp ± 0.007th U, which after evo- 
lution to the Z-boson mass scale results in 
a s (M z ) = 0.1212(11). The dominant quantifi- 
able theory uncertainty resides in as yet uncalcu- 
lated higher-order QCD corrections and improve- 
ments of the perturbative series through renor- 
malisation group (RG) methods. 

Most suitable for the a s determination is the 
t decay rate into light u and d quarks R t ,v/a 
via a vector or axialvector current, since in this 
case power corrections are especially suppressed. 



Theoretically, R r yiA takes the form lj 
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where Sew = 1.0198(6) [5] and S' EW = 0.0010(10) 
[B] are electroweak corrections, 5^ comprises the 
perturbative QCD correction, and the 5^)y, A de- 
note quark mass and higher D-dimcnsional oper- 
ator corrections which arise in the framework of 
the operator product expansion (OPE). 

When computing the total r hadronic width a 
phase-space integral over the physical spectrum 
has to be calculated, which by analyticity can be 
related to a contour integral over QCD correlators 
in the complex s-plane, where s is the invariant 
mass of the final state hadronic system. As we 
also intend to RG improve the perturbative se- 
ries, the questions arises if the RG resummation 
should be performed before or after evaluating 
the contour integral? If the "true" all-order re- 
sult were available, both treatments should agree, 
but to any finite order in perturbation theory, sig- 
nificant differences may arise. 

The approach of first RG-improving the cor- 
relators and then performing the contour inte- 
gral was introduced in [718] and termed contour- 
improved perturbation theory (CIPT), as unar- 
guably large running effects of a s (^/s) along the 
contour are resummed. However, it is known that 
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the QCD series is divergent, being asymptotic at 
best. Hence, also the explicit expansion coeffi- 
cients of the perturbative series are bound to be- 
come large. If cancellations between explicit coef- 
ficients and running effects occur, just resumming 
the running effects may not lead to a good ap- 
proximation, but it could be better to perform a 
consistent expansion in powers of a s (M T ), being 
called fixed-order perturbation theory (FOPT). 

While it can be demonstrated that, as ex- 
pected, CIPT is a good approximation, and CIPT 
as well as FOPT are compatible, when the run- 
ning effects dominate [9], in the large-/?o approxi- 
mation, in which 6^ is exactly calculable, FOPT 
gives a better approximation to the true result 
and CIPT is not seen to be compatible with it 
|10|llj . Thus it is not a priori clear which be- 
haviour prevails in real QCD. It should be obvious 
that this is a question about perturbative orders 
beyond the presently known O(a^) term [12] . 

To further investigate the difference between 
CIPT and FOPT, the contribution of higher per- 
turbative orders has to be modelled. This then 
allows to investigate the following questions: 

i) Are FO and CI perturbation theory found 
to be compatible, once terms beyond the 
currently known perturbative coefficients of 
the series for <5' ' are included? 

ii) How do FO and CI perturbation theory at a 
particular order compare to the true result 
for fll°>, and is the closest approach to the 
true result related to the minimal terms in 
the respective series? 

iii) And finally, which of the two methods, 
FOPT or CIPT provides the closer ap- 
proach to the true value at order 0{a 4 s ), 
and in general? 

The working assumption in all this is that the 
"true" result is approximated with reasonable ac- 
curacy by the Borel sum of the model series un- 
der consideration, since the power corrections to 
R T are known to be small. The model employed 
below will be constructed such as to incorporate 
the general structure of the Adlcr function in the 
Borel plane, dictated by the OPE and the RG. 



2. PERTURBATIVE CORRECTION 

In the following, we shall only be concerned 
with the purely perturbative correction which 
gives the dominant contribution to R T y/A- In 
FOPT it takes the general form 

oo n 

4o = E °( M ')" E k c ^ k Jk -i > ( 3 ) 

n=l k=l 

where a([i 2 ) = a M = a s (/i)/7r, and c 7li k are the 
coefficients which appear in the perturbative ex- 
pansion of the vector correlation function, 

M*) = -^E<E^ ln * ■ w 

n=0 k=0 V /* / 

At each perturbative order, the coefficients c n> \ 
can be considered independent, while all other 
c n ,k with k > 2 are calculable from the RG equa- 
tion. Further details can for example be found 
in ref. (T3]- Finally, the J; are contour integrals 
which are defined by 

-L I— (l-z) 3 (l + z)W(-z). (5) 
2m J x 

\x\ = l 

The first three which are required up to 0(a^) 
take the numerical values 

7—1 T — 19 T — 265 1 _2 (f-\ 

At order a™ FOPT contains unsummed loga- 
rithms of order ln'(— x) ~ n l with I < n related 
to the contour integrals J; . CIPT sums these log- 
arithms, which yields 

oo 

= £cm-Wt 2 ) (7) 

71=1 

in terms of the contour integrals J^(M^) over the 
running coupling, defined as: 

,r n {Ml) = l-j>^{l- x f{l +x )a n \-M 2 T x). 
1*1=1 

(8) 

In contrast to FOPT, for CIPT each order n 
just depends on the corresponding coefficient c n< i. 
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Thus, all contributions proportional to the coeffi- 
cient c„ i which in FOPT appear at all perturba- 
tive orders equal or greater than n are resummed 
into a single term. 

Numerically, the two approaches lead to signif- 
icant differences. Using a s (M T ) = 0.34 in eqs. ([3]) 
and (0, one finds 

4o = 0.2200 (0.2288), (9) 

4°i = 0.1984 (0.2021), (10) 

where the first number in both cases employs the 
known coefficients up to 0(a^) [T2] and the num- 
bers in brackets include an estimate of the 0{a 5 8 ) 
term with €5,1 ~ 283 [13) . Inspecting the individ- 
ual contributions from each order, up to 0{a 5 8 ) 
the CIPT series appears to be better convergent. 
However, around the seventh order, the contour 
integrals J^(M^) change sign and thus at this or- 
der the contributions are bound to become small. 
Therefore, the faster approach to the minimal 
term does not necessarily imply that CIPT gives 
the closer approach to the true result for the re- 
summed series. 

3. A PHYSICAL MODEL 

To clarify whether FOPT or CIPT results in a 
better approximation to S^°\ one needs to con- 
struct a physically motivated model for its series. 
The corresponding model will be based on the 
Borel transform of the Adler function Dy(s): 

d v (s) = -s±n v ( s ) = ^ [1+S00] . (ii) 

In the following discussion it is slightly more con- 
venient to utilise the related function D(s). Its 
Borel transform B[D](t) is defined by the relation 

00 

D(a) = J ' dte- t/a B[D](t). (12) 


The integral D(a), if it exists, gives the Borel sum 
of the original divergent series. It was found that 
the Borel-transformed Adler function B[D](t) ob- 
tains infrared (IR) and ultraviolet (UV) renor- 
malon poles at positive and negative integer val- 
ues of the variable u = 9t/(in), respectively 
|14ll5j . (With the exception of u = 1.) 



Apart from very low orders, where a dominance 
of renormalon poles close to u = has not yet set 
in, intermediate orders should be dominated by 
the leading IR renormalon poles, while the lead- 
ing UV renormalon, being closest to u = 0, dic- 
tates the large-order behaviour of the perturba- 
tive expansion. Assuming that only the first two 
orders are not yet dominated by the lowest IR 
renormalons, one is led to the ansatz 

B[D](u) = B0Y v }{u) + B[D™}{u) + 

B[D^}(u) + 4 + 4°u, (13) 

which includes one UV renormalon at u = — 1 , the 
two leading IR renormalons at u = 2 and u = 3, 
as well as polynomial terms for the two lowest per- 
turbative orders. Explicit expressions for the UV 
and IR renormalon pole terms B[Dp V ](u) and 

i?[-Dp R ](w) can be found in section 5 of ref. [13]. 

Apart from the residues v and d 1 ^, the full 
structure of the renormalon pole terms is dic- 
tated by the OPE and the RG. Therefore, the 
model (|13p depends on five parameters, the three 
residua d^ Y , d!> R an d <^3 R , as well as the two poly- 
nomial parameters dy° and d^° . These parame- 
ters can be fixed by matching to the perturbative 
expansion of D(s) up to 0(a 5 s ). Thereby one also 
makes use of the estimate for c§^\. The parame- 
ters of the model (fl~3|) are then found to be: 

d l l Y = -1.56-10~ 2 , df = 3.16, 4 R = -13.5, 

d%° = 0.781 , d\° = 7.66 • 10~ 3 . (14) 

The fact that the parameter d^° turns out to be 
small implies that the coefficient C2,i is already 
reasonably well described by the renormalon pole 
contribution, although it was not used to fix the 
residua. Therefore, one could set df° = and 
actually work with a model which only has four 
parameters. The predicted value 05,1 = 280 in 
this model turns out very close to the estimate, 
which can be viewed as one test of the stability 
of the model. 

A graphical account of the model (fT3|) for the 
(reduced) Adler function D(M^) is displayed in 
figure [T] The full circles denote the partial sums 
of the perturbative series up to order n. The 
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Perturbative order n 



Figure 1. Results for D(M^) (full circles) at 
a s {M T ) — 0.34, employing the model (fl3|) . as a 
function of the order n up to which the terms in 
the perturbative series have been summed. The 
straight line represents the result for the Borel 
sum of the series. 



minimal term of the series at the 5th order is 
marked by a grey diamond. The perturbative 
results are compared with the Borel sum of the 
model (straight line)!]] Generally, figure Q] shows 
that the model is well-behaved: the series goes 
through a number of small terms such that the 
truncated series nicely agrees with its Borel sum, 
before the sign-alternating asymptotic behaviour 
takes over around n = 10. 

The implications of the model (fT3"|) for 5^ in 
FOPT and CIPT is graphically represented in fig- 
ure [21 The full circles denote the result for S^q 
and the grey circles the one for Sqj, as a func- 
tion of the order n up to which the perturbative 
series has been summed. The straight line cor- 
responds to the principal value Borel sum of the 
series, 5gg = 0.2371, and the shaded band pro- 
vides an error estimate based on the imaginary 
part divided by tt. The order at which the se- 
ries have their smallest terms is indicated by the 
grey diamonds. As is obvious from figure [21 like 
for the Adler function itself, FOPT displays the 
behaviour expected from an asymptotic series: 
the terms decrease up to a certain order around 
which the closest approach to the resummed re- 

1 Also shown as the shaded band is an estimate of the un- 
certainty inferred from the complex ambiguity which arises 
while defining the Borel integral over the IR renormalon 
poles. For details see appendix A of ref. 1 131 . 
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Fi gure 2. Results for S F q (full circles) and <5qj 
(grey circles) at a s (M T ) = 0.34, employing the 
model (|13|) . as a function of the order n up to 
which the terms in the perturbative series have 
been summed. The straight line represents the 
result for the Borel sum of the series. 



suit is found, and for even higher orders, the di- 
vergent large-order behaviour of the series sets in. 
For CIPT, on the other hand, the asymptotic be- 
haviour sets in earlier, and the series is never able 
to come close to the Borel sum. 

The finding that in the model JT^ CIPT misses 
the full Borel sum can be traced back to the fact 
that in CIPT the running effects along the com- 
plex contour are resummed to all orders, while 
explicit contributions of the c Mi i at a certain or- 
der are dropped. However, being an asymptotic 
series, also the Adler function coefficients c n ,i be- 
come large, and cancellations between the explicit 
contributions and the running effects take place. 
As was shown in section 4 of [13], in the large- /3q 
approximation, for R T the leading IR. renormalon 
cancels completely, and also the large-order di- 
vergence of the series is softened. Even though in 
real QCD the leading IR does not anymore can- 
cel completely, for R T it is still suppressed by a 
factor 1/n 2 , and furthermore a sign-alternating 
UV renormalon component does not yet show up 
in the known coefficients. Thus, the cancellations 
between running effects and explicit coefficients 
are also expected to prevail in full QCD. 

The deficiency of CIPT to approach the Borel 
sum of the series, which leads to the marked dif- 
ferences of (5pQ and 5^, can also be observed 
when the Adler function is inspected along the 
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complex circle. This is discussed in detail in ap- 
pendix B of ref. [13] ■ While FOPT converges to 
the Borel sum on the full circle (though rather 
badly close to the Minkowskian axis), in some re- 
gions of the circle CIPT largely differs from the 
resummed result, again due to the missed cancel- 
lations. 

As the behaviour of CIPT versus FOPT hinges 
on the contribution of the leading IR renormalon 
at u = 2, in principal also models can be con- 
structed for which CIPT provides a good ac- 
count of the Borel sum. These would generally 
be models where d!> R 1S much smaller than the 



value quoted in eq. ([14]) . While such models can 
at present not be excluded, the pattern of the 
individual contributions appears more unnatural 
than in the model (|13| : the known c n ^\ can only 
be reproduced when one allows for large cancel- 
lations between the individual terms. Thus, the 
behaviour generally expected from the presence 
of renormalon poles, namely dominance of lead- 
ing IR poles at intermediate orders, would be lost. 

4. DETERMINATION OF a s 

The starting point for a determination of a s 
from hadronic r decays is eq. © for the decay 
rate of the r lepton into light u and d quarks. The 
analysis will be based on FOPT together with 
the ansatz (|T3| for higher-order terms discussed 
in the last section. Due to the observation that 
CIPT is not able to approach the resummed se- 
ries, it will not be employed below. 

The first step of the a s analysis consists in 
estimating the values of the power corrections 

^ud V+A' which arise from higher-dimensional op- 
erators in the framework of the OPE. Given 
these estimates and experimental data, a phe- 
nomenological value of 5^ can be calculated us- 
ing eq. This allows to determine the value 
of a s (M T ) by requiring that the theoretical value 
(5pQ matches the phenomenological value S^ cn . 
Errors are estimated by varying all parameters 
within their uncertainties. 

In view of the smallness of the light quark 
masses m u and nid, as well as the suppression of 
the dimension-4 contribution in R T , the dominant 
power correction arises from the six dimensional 



4-quark condensates. As the number of contribut- 
ing operators is too large to treat the 4-quark con- 
densates individually, conventionally the so-called 
vacuum-saturation approximation (VSA) [16] is 
employed. Then the corresponding contribution 
takes the form 



°(qqqq),V+A 



512 



n 3 a. 



27 " "* M$ 
(-4.8 ±2.9) • 1CT 3 . 



(15) 



where for the numerical estimate the required 
quark condensate is taken from the GMOR re- 
lation [17118] , and p = 2 ± 1 was assumed to take 
into account violations of the VSA. This choice 
includes most estimates of the four-quark conden- 
sate present in the literature. 

To complete the estimate of power corrections 
to R T , the longitudinal contribution which arises 
from the pseudoscalar correlator still has to be 
includcdo Because the perturbative series for 
this correlator does not converge very well, the 
approach of refs. |19|20j will be followed. The 
main idea is to replace the QCD expression for 
the pseudoscalar correlator by a phenomenologi- 
cal representation. The dominant contribution to 
the pseudoscalar spectral function stems from the 
well-known pion pole, giving 



C7T 

°ud,S+P 



= -16tt 



f 2 M 2 

Mi 



Ml 

Ml 



(16) 



plus small corrections from higher-excited pionic 
resonances. Repeating the analysis of ref. [19] and 
updating the input parameters, one finds 



Jud,S+P 



(-2.64 ±0.05) ■ 10" 



(17) 



Collecting all contributions, and adding the errors 
in quadrature, one arrives at the total estimate of 
all power corrections: 



^PC 



(-7.1 ±3.1) -10" 



(18) 



The value {T5J is consistent with the most recent 
fit to the r spectral functions performed in ref. [4] . 

As a matter of principle, the OPE of correla- 
tion functions in the complex s-plane could be 



2 The scalar correlator, being proportional to (m u — m d ) 2 , 
is completely negligible. 
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inflicted with so-called "duality violations" [2"T] . 
These arise from the contour integral close to the 
physical region which even though suppressed in 
R T could be sizeable [52]. Nevertheless, before a 
possible additional duality violating contribution 
can be extracted consistently from a combined fit 
to spectral moments, it shall be omitted. 

Employing the value R T y+A = 3.479 ± 0.011, 
which results from eq. ([T]) in conjunction with 
R r S = 0.1615 ± 0.0040 g], as well as \V ud \ = 

0. 97418 ±0.00026 [23], from eqs. © and (HHJ) the 
phenomenological value for 5^ can be derived: 

5^ en = 0.2042 ± 0.0050 . (19) 

The dominant experimental uncertainty in (I19p 
is due to R T y + A and the theoretical one to the 
dimension-6 condensate. The final step in the 
extraction of a s (M T ) now consists in finding the 
values of a s for which S^ cn matches the theoret- 
ical prediction, which yields [13] 

a s (M T ) = 0.3156±0.0030 cxp ±0.0051 th . (20) 

Evolving this result to the Z-boson mass scale 
finally leads to 

a s (M z ) = 0.11795 ±0.00076, (21) 

in perfect agreement with the world average [24] . 
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